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Abstract 

The dynamics of two coupled spins- 1/2 coupled to a spin-bath is studied as an extended model of 
the Tessieri-Wilkie Hamiltonian pj. The pair of spins served as an open subsystem were prepared 
in one of the Bell states and the bath consisted of some spins- 1/2 is in a thermal equilibrium state 
from the very beginning. It is found that with the increasing the coupling strength of the bath 
spins, the bath forms a resonant antiferromagnetic order. The polarization correlation between 
the two spins of the subsystem and the concurrence are recovered in some extent to the isolated 
subsystem. This suppression of the subsystem decoherence may be used to control the quantum 
devices in practical applications. 
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I. INTRODUCTION 



Quantum decoherence is a common and inevitable phenomenon in an open quantum 
system due to its interaction with the surrounding bath or environment. An initial 
superposition state of the system, ps{t = 0) = \ipo){'^o\, has to decay into a classical mixture 
of states, ps = J2iPiPiy J2iPi = 1> after a decoherence time r^. During the same period 



of time, some of t 



he m: 



'ormation in the initial state of the system might irreversibly lose 



into the bath [2], [s, Q, [s, 0]. A real-world system, for instance a quantum device or qubit, 
cannot be completely isolated from the environment. In the last two decades, there are 
great interests in the search for realizations of quantum computation as well as quantum 
communications, such technologies rely on the possibilities that quantum devices can be 
manufactured with negligible decoherence in the period of operation. Since decoherence 
is intrinsic to open quantum systems, the problem is transformed to a problem on how to 
reduce or eliminate the decoherence of an open quantum system. Many works have been 
devoted into the research of the influences caused b^ the subsystem-bath coupling and 
the decoherence dynamics of open quantum systems 
quantum communication and quantum computation. 



la la 



for the possible realizations of 



A very important part of any theoretical research of decoherence is the modeling of 
the bath or environment. There had been mainly two important types of bath schemes: 
(i) the bath consisted of harmonic oscillators, e.g., the Caldeira-Leggett model Js]; (ii) 
the bath made up of spi.s-1/2. e.g.. model used in Ref. M. For both types of bath 
model, there are typically two kinds of approximations, Markovian or non-Markovian 
[1^ . used in the studies. The focus of most of the researches were on the eliminating 
the destructive effects of the environment to the system. However, Tessieri and Wilkie 
proposed a new idea to the reducing of decoherence of the open subsystem, which is 
a spin-1/2 coupled to a low-temperature bath of spins. They introduced the coupling 
terms between bath spins into their spin-bath Hamiltonian, which is a simplification of 
the vibronic Hamiltonian of the impurity and crystal [l^ using spin-1/2 modes. One of 
the most important results they found is that the strong antiferromagnetic intra-bath 
interactions could be utilized to make the dynamics of the central spin almost au- 
tonomous from the bath around. And Dawson et al. 141 also found the decoherence effect 
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of the bath can be suppressed by the increasing intra-bath couphng in Tessieri-Wilkie model. 



Inspired by their works and the papers of Xiang et al. 



and Liu et al. jlj!] , the aim of 



this paper is to study the dynamics of correlation characteristics between the two coupled 



center spins, w^ 
by the Ref. [l 



lich constitute our subsystem. The bath in our study is the same model used 
[l^ . which is prepared in a thermal equilibrium state at time t = 0. The 



dynamics of the subsystem and the bath then calculated, phenomena such as decoherence 
oscillation were observed, which can be quantified by the amplitudes of polarization com- 
ponents {(Jx), {(Jy) and {(Tz)- We studied the evolution of the spatial correlation between 
polarizations of the two subsystem spins along three directions x, y and z. we also dis- 



171] , which is an essential ingredient 



cussed the entanglement between the two system spins 

in the quantum communication jisl . and has no classical counterpart. The concurrence 
of the subsystem was chosen to be a measure of the entanglement between the two spins 



within the subsystem [20|, 1211]. It will be demonstrated that in our extended model, both 
the polarization correlation and entanglement of the subsystem can be restored to a great 
extent to the isolated case by the increasing pairwise couplings between the bath spins. And 
we also try to clarify the physics behind this revival function of intra-bath coupling. The 
rest of this paper is organized as follows. In section [Tll we introduce the Hamiltonian for 
our two-center-spins-spin-bath model; In section UTTt we explain the computation procedures 
about the evolution of the correlation and of the concurrence for the subsystem; Detailed 
results and discussions are in section [IVl The conclusion of our study is given in section |Vl 



II. THE EXTENDED SPIN-SPIN-BATH MODEL 



The subsystem we studied here consists of two spins antiferromagneticly coupled in 
the X direction, and aligned in x and z direction respectively by applied fields. The bath 
consists of — 2 spins, every pair of the spins are also coupled antiferromagneticly in 
the X direction, and aligned in both x and z direction. The subsystem interacts with the 
bath by the antiferromagnetic coupling in the x direction. This model is an extension to 
the Tessieri-Wilkie model |l]], which comes from the simplification of an experimentally 
realizable model system. The original model is based on the system that an interstitial He 
atom in an otherwise pure silicon or diamond cluster, the He atom causing small lattice 
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distortion in the cluster 



22], thus inducing vibronic couphng between the electronic states 



and phonon states of the cluster. This mechanism is further simplified by considering only 
the ground state and the first excited state of the interstitial He atom, representing the 
states with a spin- 1/2 Pauli matrix, and representing the phonon states as spin-1/2 modes, 
which resulting in the Tessieri-Wilkie model. This model system is regarded as a potential 
realization for quantum computing device. Our extension is that instead of one center 
spin, we have two center spins coupled together by an antiferromagnetic coupling. In this 
extension we merely regard the Tessieri-Wilkie model as a model system and did not go 
into details of the physical realization. So that we will only focus our attention to the 
extended model itself in the following discussions. 



The Hamiltonian of our system can be written as: 

H = Hs + Hb + Hsb, (1) 
Hs = f crW + + + ^a^ + X^^a^^^a^^^ , (2) 

N N N-l N 

i=3 i=3 1=3 j=i+l 

N 

Hsb = A.,5^(cr«a» + ai^)a«). (4) 

j=3 

Where Hs is the Hamiltonian of the subsystem with two interacting spins (labeled by 1,2) 
coupled in x direction with coupling strength A^^. /3, and Uh may be regarded as applied 
fields acting in x and z directions to the sub-system and bath spins. The ax and are the 
Pauli matrices: 









fl \ 
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(5) 



The bath consists of — 2 spins labeled by 3 to N, coupled each other in the x direction 
with coupling strength Af,?,, denoted as Hb- The coupling between subsystem and bath is 
given by Hsb, which is also the coupling of each subsystem spin with every bath spin in x 
direction, with coupling strength As^. All the interactions are antiferromagnetic, that is 

A.s>0, Xbb>0, Xsb>0. (6) 

In the following calculation, we use the system of units where the Boltzmann constant 
kB = ^ and Plank constant h = 1. The other parameters in the Hamiltonian take the 



following values in the subsequent computations: 

uj^ = 0.8, (3 = 0.1, Uh = 1.0. 



III. CALCULATION PROCEDURES 

In this section, we describe in detail the algorithms for the calculation of the dynamics of 
the subsystem modeled by the Hamiltonian ([T]). The initial state of the subsystem is chosen 
to be one of the four Bell states: 

1^5(0))^ = l/v^(|ll) + |00)), (7) 
1^5(0))^ = 1/V2(|10) + |01)), (8) 
1^^5(0))^ = 1/^2(111) -|00)), (9) 

IM^r = i/v^(|io)-|oi)). (10) 

where |0) and |1) refer to the spin "down" and spin "up" in the z direction, respectively. 
The reason that we use these states as initial states is that among all the quantum states 



for a pair of coupled spin- 1/2, Bell 
other interesting characteristics 



_ Dases have the largest degree of entanglement and many 
231]. However, it should be noted that the last Bell state 
given in equation (|TOl) is a state free from the bath, the action of the interaction part of the 
Hamiltonian, Hsb, on the state gives zero. So that the state will evolve with time just as an 
isolated system. Further more, the state is also an eigenstate of the subsystem Hamiltonian 

H,{l/V2{\10) - |01))) = -l/v^A,,(|10) - |01)). (11) 

The time evolution of the state is simply an exponential factor e*"^*, with u = —Xss- So that 
the physical quantities will not vary with time in this state and we will not consider it in the 
following calculations. The initial state of the bath is taken to be the thermal equilibrium 
state: p_b(0) = exp(— if^/T), where Z is the partition function of the environment 
Z = Tr{exp{—HB/kBT)). The density matrix p{t) of the whole system can formally be 
represented as: 

p{t) = ex-p{-iHt)p{0)exp{iHt) (12) 
p(0) = P5(0)®Pb(0) (13) 
P5(0) = \MO)){MO)\ (14) 



n order to find the density matrix p{t), we follow the method suggested by Tessieri et al 
l|. The thermal bath state Pb(0) can be expanded with the eigenstates of the environment 
Hamiltonian: 
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Pb{0) = 



n=l 



n=l 



-Er,/T 



(15) 

(16) 
(17) 



Here |0n^^), n = 1, 2, 3, ■ ■ ■ , 2™, m = N — 2, are the eigenstates of the environment Hamilto- 
nian Hb, and En the corresponding eigen energies. With this expansion, the density matrix 
p{t) can be written as: 

2™ 

p{t) = J2^n\'i^nit)){^n{t)\. (18) 
n=l 

Where 

|^„(t)) = exp(-2m)|^„(0)) = [/(t)|^„(0)). (19) 

The initial state is 

|vI/„(O)) = |V^5(O))|0f)). 
The evolution o per ator U(t) can be evaluated by the efficient algorithm of polynomial 



schemes 



24 



25 



26|. The method used in this calculation is the Laguerre polynomial ex- 



pansion method we proposed in Ref. [26|, which is pretty well suited to this problem and 
can give accurate result in a comparatively smaller computation load. More precisely, the 
evolution operator U (t) is expanded in terms of the Laguerre polynomial of the Hamiltonian 



as: 



U{t) 



-iHt 



a+1 oo 



it 



Ll{H) 



l + it ^ \ l + it 

where a distinguishes different types of Laguerre polynomials j27|, k is the order of the 
Laguerre polynomial. In real calculations the expansion has to be cut at some value of 
kmax, which was taken to be 24 in this study. With the largest order of the expansion fixed, 
the time step t is restricted to some value in order to get accurate results of the evolution 
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TABLE I: the number of the bath states, M, needed in the calculation for different temperatures 
and thresholds of uJth- 



Temperature T 


< 0.04 


0.05 


0.08 


0.10 


0.15 


0.20 


M{uJn > 0.00001) 


8 


8 


28 


70 


70 


70 


Micjn > 0.0001) 


1 


8 


8 


28 


28 


70 


M{uJn > 0.001) 


1 


1 


8 


8 


28 


70 



operator. For longer times the evolution can be achieved by more steps. The action of the 
Laguerre polynomial of Hamiltonian to the states is calculated by recurrence relations of 



26| is about 8 times as 



the Laguerre polynomial. The efficiency of this polynomial scheme 
that of the Runge-Kutta algorithm used in Ref. When the states \'^n{t)) are obtained, 
the density matrix can be obtained by performing a summation in equation f|T8l) . 

When the temperature is not very high, the weights for high energy states will be very 
small so that only few lowest states need to be considered. The maximum states to be 
included in the summation can be determined by setting a threshold of weight uth, and 
keep only those states whose weight Un > ujth- As an example, table ([T]) gives the number 
of states to be used for different temperatures and threshold. In this case the number of 
environment spins m = 8, and couplings among the environment spins Xbf, = 4.0. As can 
be seen from the table, the number of states needed is much less than the total number of 
states, 2® = 256. In real calculation the up bound of the summation in equation (fTSll will 
be M << 2™, chosen by the criteria specified for the accuracy of results. That is to say, 
equation (|T8l) could be changed into the following equation: 

M 

p(t) = ^cu„|^„(t))(*„(t)|. (20) 

n=l 

After obtaining the density matrix of the whole system, the reduced density matrix is 
calculated by a partial trace operation to p(t), which trace out the degrees of freedom of the 
environment: 

psit)=TrBipit)). (21) 

For the model of this paper, ps{t) is the density matrix of the open subsystem consists 
of two center spins, which can be expressed as a 4 x 4 matrix in the Hilbert space of the 



subsystem spanned by the orthonormal vectors |00), |01), |10) and |11). 



After the reduced density matrix is obtained, any physical quantities of the subsystem 
can be obtained easily. In the following we concentrate on two important physical quantities 
of the subsystem which reflect the decoherence and entanglement degree of the subsystem 
state. These two quantities are the spatial polarization correlation and the concurrence. 



A. Polarization correlation 



The polarization of either spin in the subsystem is defined as 

P«(t) = Tr(p«(t)a«), i = 1, 2, (22) 

where a — axi + (Jyj + azk. To simplify the calculation, we may also trace out one of the 
spin degrees of freedom to obtain the 2x2 sub-reduced density matrix for each spin: 



p«(i)oi p«(^)oo 



Here Trj means to trace out the other degrees of freedom of i, i.e. trace out 2 when i — 1 
and trace out 1 when i — 2. Then the three components of P^'^\t) are expressed as: 

«.) = Tr (p« W^x.) = P«(i)io + P«(t)oi, (23) 

{al) = Tr (p^HtK) = ^ {p^HtU - P^Ht)oi) , (24) 

{ai) = Tr (p«(^)^.) = P^Ht)n - P«(^)oo. (25) 
The polarization P*^') (t) may be viewed as an indicator of quantum decoherence. 



The correlation between the two subsystem spins can be described by the correlation 
functions defined bellow: 

C.. - (26) 

Cyy - (27) 

C.. = (28) 
where (aa^aa^) = Tr (^ps{t)aa^ aa^^ , ck — x,y,z. Since — ay = = 1, so that 

^(Ta ^ + aa^^ — 2 + 2aa^aa \ the correlation function is thus also a measurement of the 
fluctuations of the total spin of the system. 
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B. Concurrence 



The concurrence of the two spin-1/2 system is an indicator of their intra entanglement, 



which is defined as 



20|: 



C = max{Ai - A2 - A3 - A4, 0}, (29) 

where Aj are the square roots of the eigenvalues of the product matrix psPs in decreasing 
order. Equation fl2^ applies to all kinds of states, either mixed or pure. The matrix is 
constructed as (cTj, (g) o-y)p*s{o'y ® o-y)- If the bipartite quantum state ps is pure |28|, such as 
the states in equations (17|)- (fT0|l . They can be written as: 



Ps = 

1^) = a|00) + 6|01) +c|10) + 



then equation fl29l) could be simplified to 



C(|^)) = 2\ad-bc\. 



(30) 



IV. RESULTS AND DISCUSSIONS 



In this section we give the calculated results of polarization correlations and concurrence 
with discussions. In all the calculations given here, the parameters A^;, = A^s = 1-0 unless 
otherwise specified. And the temperature is set as T = 0.1 much higher than T = 0.02 in 
Ref. pj] in order to stress the generality of the algorithm and conclusion. The results for 
isolated systems Xgb = 0.0 are presented as a standard for comparison, in which there is 
no decoherence occurs. On the other hand, the case of Xbb = 0.0 will also be considered, 
in which there is a strong decoherence occurs and the subsystem initial state is beyond 
retrieval. 



A. Polarization correlation 

The polarization correlations of the isolated system are plotted in figures [H [2] and [31 we 
see that the evolutions of the polarization correlations are periodic in time. Then we take 
into account the case that Xbb = 0.0, which means there is no coupling among bath spins, 
in order to highlight the effect of strong intra-bath coupling. In figure HJ we show the pure 
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time 



FIG. 1: Evolution for polarization correlation of the isolated subsystem, Caa = (ca ctq ) — 
a is X, y or 2; , between the two spins in the subsystem. The initial state of the 
subsystem is l/x/2(|00) + 




FIG. 2: Evolution for polarization correlation of the isolated subsystem, C^a = {o'a^f^a^) — 
a is X, y or z, between the two spins in the subsystem. The initial state of the 
subsystem is 1/V2(|01) + |10)). 
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time 

FIG. 3: Evolution for polarization correlation of the isolated subsystem, Caa = 

a is X, y or z, between the two spins in the subsystem. The initial state of the 
subsystem is l/^/2{\m) - |11)). 

destruction effect on polarization correlation of subsystem spins imposed by the bath spins 
(m = 6) to the open subsystem, where we neglect the intra-bath coupling strength. The 
initial state of the subsystem is |V^s(0))^ = 1/72(101) + |10)). 



Figure [5] through figure [13] are the polarization correlations when the subsystem coupled 
to the bath with different X^b- The corresponding results for isolated systems are also drawn 
on each figure as solid lines. In these figures, we have m = 4 spins in the bath, though the 
number of bath spins is small, it still can give some information on the influence of bath cou- 
pling to the system. We also did calculations with more bath spins, m = 6, show in figure [TH 



It is clear from the figures that when the subsystem coupled to the bath, all the polariza- 
tion correlations are different from the isolated systems. The difference reflects the influences 
of the bath and decoherence of the subsystem. However, as pointed out in reference [l| and 



14| for the case of a single spin system, the coupling strength among environment spins can 



suppress the decoherence and partially restore the subsystem coherence. Figure [5] through 
figure [13] are plots of the polarization correlations as function of the time with different 
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5 10 15 20 25 30 35 40 

time 



5 10 15 20 25 30 35 40 

time 



(a)(a>2) _ A.^ = 0.0 (h){alal) ~ {al){al), A,, = 0.0 




5 10 15 20 25 30 35 40 



time 

(c)(aia2)-(ai)(a2), Ab6 = 0.0 

FIG. 4: The evolution of polarization correlation along three directions of the open subsystem. 
Where there are 6 spins in the bath and the initial state of the subsystem is l/-v/2(|10) + 101)) and 
Ass = Asfe = 1.0, Xbb = 0.0. 

initial states and increasing intra-bath coupling strength. At weak intra-bath coupling, the 
results are drastically different from the isolated systems and similar with figure HI as the 
intra-bath coupling is increased, the difference between the open subsystem and the isolated 
system becomes smaller; and at very strong strengths, the difference becomes very small 
especially in the case of \ips{0))^. The other two initial Bell states will be dissipated when 
time is elongated. 

These results suggest that by changing the bath coupling strength, one can effectively 
control the subsystem coherence properties. The suppression of decoherence effect is also 
depends on the number of bath spins, it is found that the more bath spins, the less effective 
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TO 




- >.^^=1 .0 

. X. , =1.0. X..=2.Q 




ml 



- k^=^ .0 

. k , =1.0, X..=4.0 



(a)Ahb = 2.0 





- '.33=1.0 



{c)Xbb = 6.0 



(b)A6b = 4.0 




. k ^=1.0,?, =8.0 



20 25 

time 



30 35 40 



(d)Ab6 = 8.0 



FIG. 5: Evolution for polarization correlation along x direction of the open subsystem, {cr^cr'^) — 
{^x){^x)- The initial state of the subsystem is l/\/2(|00) + |11)). 

of suppression of the decoherence with the same bath couphng strength. To show this, we 
provide the results of Xbb = 6.0 and Xbb = 10.0 on the Bell state of lipsiO))"^ = 1/V2(|01) + 
1 10)) in figure [m where there are m = 6 bath spins. The values of other parameters are 



the same as before. Through the comparison of figure 14(a) with figure 8(c), we see that 



when m = 4, Xbb = 6.0 is almost sufficient to recover the isolated case except a slightly 



elongated period (see figure 8(c)), but when m is increased to 6, the coupling strength 
need to be increased at least to Xbb = 10.0 to get the same effect (see figure 14(b)). The 



same feature also happens to the evolution of Cyy and C^z which can be seen by comparing 



figure 14(c) with figure 9(c) and figure 14(e) with figure 10(c) If we go on to increase the 
coupling strength Xbb, we can suppress the decoherence and disentanglement effect due to 
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5 10 15 20 25 30 35 40 5 10 15 20 25 30 35 40 

time time 



(c)Af,b = 6.0 {d)Xbb = 8.0 

FIG. 6: Evolution for polarization correlation along y direction of the open subsystem, (ayay) — 
{al){al). The initial state of the subsystem is l/^/2(|00) + |11)). 

the bath with more spins. In figure [T5|, we compare the results of the cases of \hb = 4.0 and 
Abb = 24.0, in which the subsystem seems to be decoupled from the bath; the initial state is 
also l/v/2(|01) + |10)). 
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5 10 15 20 25 30 35 40 5 10 15 20 25 30 35 40 

time time 



(c)A66 = 6.0 (d)A66 = 8.0 

FIG. 7: Evolution for polarization correlation along z direction of the open subsystem, (cr^a^) — 
(cr^)((T2). The initial state of the subsystem is 1/V2(|00) + 
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time time 

{c)Xbb = 6.0 (d)A66 = 8.0 

FIG. 8: Evolution for polarization correlation along x direction of the open subsystem, (cr^cr^) — 
{al){al). The initial state of the subsystem is 1/V2(|01) + |10)). 
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(a)A, 



66 



2.0 



(b)A 



66 



4.0 




time 




{c)Xbb = 6.0 



{d)Xbb = 8.0 



FIG. 9: Evolution for polarization correlation along y direction of the open subsystem, {CyTy) 
{al){al). The initial state of the subsystem is 1/V2(|01) + |10)). 
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(a)A, 



66 



2.0 



{h)Xbb = 4.0 




{c)Xbb = 6.0 



{d)Xbb = 8.0 



FIG. 10: Evolution for polarization correlation along z direction of the open subsystem, (cr^cr^) 
{al){a'^). The initial state of the subsystem is 1/V2(|01) + |10)). 
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(c)A66 = 6.0 



30 35 




15 20 25 30 35 40 

time 



(d)A66 = 8.0 



FIG. 11: 



Evolution for polarization correlation along x direction of the open subsystem, (u^it^) 
The initial state of the subsystem is l/-\/2(|ll) — |00)). 
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time 



(c)A, 
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time 



(b)Afe6 = 4.0 




f 11 




- X =1.0 



15 20 

time 



(d)A66 = 8.0 



35 40 



FIG. 12: 



Evolution for polarization correlation along y direction of the open subsystem, (c^cr^) 
The initial state of the subsystem is l/\/2(|ll) — |00)). 
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10 15 20 25 30 35 40 

time 



{a)Xbb = 2.0 



. X =1 .0, X =6.0 



10 15 20 25 30 35 40 

time 



{c)Xbb = 6.0 



. X =1 .0, X =4.0 



10 15 20 

time 



{h)Xbb = 4.0 



30 35 40 



15 20 

time 



{d)Xbb = 8.0 



. X =1.0,X =8. 



30 35 40 



FIG. 13: Evolution for polarization correlation along z direction of the open subsystem, (crlaj) 
{(7l){al). The initial state of the subsystem is l/\/2(|ll) - |00)). 



21 



10 15 



30 35 



time 



{B.){alal)-{al){al),X,, = Q.Q 




10 15 



-J. =1.0 

30 35 



time 



malal)-{al){al),X,, = lQ.Q 




(d)(a>2) _ Afcfa = 10.0 




35 40 



time 



time 



(e)(aia2)-(ai)(a2), A^b = 6.0 



(f)(c7ia2)-(ai)(a2), A6fe = 10.0 



FIG. 14: Evolution for polarization correlations along three directions of the open subsystem. 
Where there are 6 spins in the bath and the initial state of the subsystem is l/\/2(|01) + |10)). 
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. X =1 .0, ?L_=4.0 



- X =1.0 



time 



(a)(^T>^>-((Ti)(a2), Ab6 = 4.0 



(b)(aia2)-(ai)(a2),A(,6 = 24.0 
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(c)Kc72)-(0(a2), A(,fa=4.0 



(d)(^>P - Afcfa = 24.0 




time 



time 



(e)(aia2)-(ai)(a2), A^b = 4.0 



(f)(c7ia2)-(ai)(a2), A6fe = 24.0 



FIG. 15: Evolution for polarization correlations along three directions of the open subsystem. 
Where there are 8 spins in the bath and the initial state of the subsystem is l/\/2(|01) + |10)). 
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B. Concurrence 




time 

FIG. 16: Evolution of the concurrence of the two subsystem spins for the isolated subsystem. The 
initial state of the subsystem is l/\/2(|00) + |11)). 

The effect of coupling strengths of bath spins on the concurrence are also studied. In 
figures [161 [13 s-^id [18] we plot the time evolution of the concurrence of three Bell states for 
the isolated subsystem. An example that the subsystem-bath coupling is considered but 
the intra-bath coupling strength is zero is showed by figure [191 which can be compared with 
figure [20] because of the same compaction condition except Xbb- 

Figures [201 [22] are plots of the evolution of concurrence of the three Bell states 

of the open subsystem. The number of bath spins is m = 6 and two coupling strengths of 
bath spins, Xbb = 6.0 and Xbb = 10.0, are used in the calculation. By comparing the results 
of Xbb = 6.0 with that of Xbb = 10.0, there are visible improvements of the concurrence by 
increasing the coupling strength of the bath spins. However, it is found that the influence 
of the same coupling strengths of bath spins has less effect on the concurrence than on the 
polarization correlations. From the figure we see that even for the strength Xbb = 10.0, the 
concurrences are still somewhat different from the isolated cases. It is also noted that the 
influence is different for different initial stats, it can be seen from figure [20] that for the 
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5 10 15 20 25 30 35 40 

time 



FIG. 17: Evolution of the concurrence of the two subsystem spins for the isolated subsystem. The 
initial state of the subsystem is l/\/2(|01) + |10)). 

case of Bell state 1^/^5(0))^ = l/v^(|01) + |10)), much better suppression of decoherence is 
observed than the other two states. Thus we provided the results for a larger bath with the 
initial state lipsiO))"^ in figure [23l where there are m = 8 bath spins. We found as Xbb = 24, 
the evolution of concurrence approaches the dynamics of subsystem evolving in isolation. 
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FIG. 18: Evolution of the concurrence of the two subsystem spins for the isolated subsystem. The 
initial state of the subsystem is 1/a/2(|00) — |11)). 




FIG. 19: The evolution of concurrence of the open subsystem. Where there are 6 spins in the bath 
and the initial state of the subsystem is l/-\/2(|10) + |01)) and Xgs = Xgb = 1.0, A;,;, = 0.0. 
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(a)A66 = 6.0 (b)A66 = 10.0 

FIG. 20: The evolution of concurrence of the open subsystem. Where there are 6 spins in the bath 
and the the initial state of the subsystem is l/-\/2(|01) + |10)) and Xgg = \sb = 1-0 
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(a)A66 = 6.0 (b)A6fe = 10.0 

FIG. 21: The evolution of concurrence of the open subsystem. Where there are 6 spins in the bath 
and the initial state of the subsystem is l/\/2(|00) + |11)) and Xgg = ^sb = 1-0 
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(a)A66 = 6.0 (b)A66 = 10.0 

FIG. 22: The evolution of concurrence of the open subsystem. Where there are 6 spins in the bath 
and the initial state of the subsystem is l/-\/2(|00) — and Xgg = Xgb = 1-0 
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(a)A66 = 4.0 (b)A6fe = 24.0 

FIG. 23: The evolution of concurrence of the open subsystem. Where there are 8 spins in the bath 
and the initial state of the subsystem is l/\/2(|01) + |10)) and Xgg = Xgb = 1-0 
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C. Discussion 



All the results of comparison in the above suggest that if we did not consider the intra- 
bath coupling, or we cannot control the intra-bath coupling, the decoherence occurring in 
our subsystem due to the bath is very severe. Although the numbers of bath spins in our 
simulation is not very big, but these bath spins can be regarded as the nearest-neighbors 
to our open subsystem, the interaction between them Xgh is much larger than the coupling 
exists between subsystem spin and other degree of freedom in the real- world environment. 
So the bath we considered in this study could form a safeguard device of the open subsystem, 
which intra-coupling could be adjusted to a high level to counteract the dissipation by itself 
and the real- world environment. 

To understand the physics behind the decoherence suppression effects due to the strong 
coupling among bath spins, we investigated the states of the bath when the coupling 
strength changed. The bath spins coupled antiferromagneticly with each other in the x 
direction, which cause the bath spins in a kind of frustrated state. The other terms in the 
Hamiltonian favors a aligned ordered state. The final state is the competition between 
the different terms as well as the thermal fluctuations and turns out to be very complex. 
However, we believe that the decoherence suppression effects are somehow related to the 
state ordering of the bath spins. 

At temperatures not high enough, the properties of the bath are determined by the few 
lowest energy levels, so we will concentrate on the lowest energy levels and try to figure out 
the characteristics of the states. In our representation where the z component of the spin 

is diagonal, it is hard to see the ordering properties of the states. Since the coupling is in 
the X direction, we thus transform our states to the x component diagonal representation 
for clarity. This is simply achieved by the following recipe to each spin state. 




1 



(31) 



Where J7 is a 2 x 2 matrix defined as 




(32) 
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Now we consider a system with 6 bath spins and two spins in the subsystem, 
the bath states can be expanded with the product state of the form \i1i2i3i4i5i6) x — 
Ih) x\'i2) xlis) x\h) x\h) xl^e) X, here ia equals or 1, the subscript x means the x-component 
diagonal representation, and a = 1, ■ ■ ■ , 6 are the index of the 6 bath spins. There are 
2^ = 64 product states, and for each product state, the subsystem can be in four states 
|00)j;, |10)^,|01)a; and \ll)x- Since the Hamiltonian Hb is invariant under permutations of 
the bath spins, so that the states with the same number of "up" spins in the x direction has 
the same expectation values of energy, thus we can group the states with the same number 
of "up" spins together. And using the number of "up" spins n to represent such states and 
denote it as \ \n). 
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gl 1 , , 1 , , , , , 1 

-1 0123456789 

FIG. 24: The probability that the bath state has half of the spins in the "up" state and half in the 
"down" state as function of the coupling strength Af,?,, the total number of bath states is m = 6. 
The initial subsystem state is l/\/2(|01) + |10)) and temperature T = 0.1. 

The probability that each kind of the product states appear in the bath state can be 
determined from the full density matrix p{t), i.e. 

P(n) = Trs{n\ \p{t) \ \n) = Tr, u;^ | (n| | vl/^(t)) |' . (33) 

The Tvs means to trace out the subsystem degrees of freedom. In the case of 6 bath spins, we 
calculated the probabihties for each state | \n) and found that the state 1 13) , which has 3 spins 
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"up" and 3 spins "down" in the x direction, has the majority probabihty, and the probabihty 
increases with the couphng strength. Figure [Ml is the plot of the probabihty -P(3). We see 
that for small coupling strength, the probability is about 0.1, and as the coupling strength 
increases, the probability increases monotonically and reaches 0.8 at A;,fe = 10, where the big 
suppression effects was obtained. Based on this observation we conclude that as the coupling 
strength increases, the bath spins are self organized to a resonant antiferromagnetic ordered 
state in x direction, one half of the spins are in the "up" state and the other half are in 
the "down" state. It is this ordering that brings the subsystem to the more coherent state. 
This mechanism may play an important role in the controlling of the subsystem coherence 
of quantum device, despite the inevitable influences of thermal noise upon the quantum 
device, we may also couple the quantum device to a system as described in this paper. By 
changing the coupling strength of the added bath spins, one may reduce the decoherence to 
a level for practical applications. 

V. CONCLUSION 

In this paper, we extended the one-center-spin-spin-bath model to a 2-center-spin- 

spin-bath model, which could be thought as an analog to two correlated qubits in quantum 
computer. By calculating the polarization correlation and concurrence of the subsystem, we 
found that the spin-bath can play a revival role in the evolution of polarization correlation 
and entanglement between two subsystem spins. In the process of calculation, we combine 
the techniques of Ref. [l| and Ref. [26| to reduce the computer resources greatly. The 
physics of this suppression was found to be the effect of the antiferromagnetic ordering of 
the bath spins in x direction. We suggest that the results may be of use in the controlling 
of decoherence of quantum devices. 

This work is supported by the National Nature Science Foundation of China under grant 
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